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ABSTRACT 

We present the results of an ongoing weak lensing survey conducted with the Subaru telescope whose initial 
goal is to locate and study the distribution of shear-selected structures or halos. Using a Suprime-cam imaging 
survey spanning 21.82 deg^, we present a catalog of 100 candidate halos located from lensing convergence maps. 
Our sample is reliably drawn from that subset of our survey area, (totaling 16.72 deg^) uncontaminated by bright 
stars and edge effects and limited at a convergence signal to noise ratio of 3.69. To validate the sample detailed 
spectroscopic measures have been made for 26 candidates using the Subaru multi-object spectrograph, FOCAS. 
All are confirmed as clusters of galaxies but two arise as the superposition of multiple clusters viewed along 
the line of sight. Including data available in the literature and an ongoing Keck spectro scopic campaign, a total 
of 41 halos now have reliable redshifts. For one of our survey fields, the XMM LSS (iPierre et aljljool field, 
we compare our lensing-selected halo catalog with its X-ray equivalent. Of 15 halos detected in the XMM-LSS 
field, 10 match with published X-ray selected clusters and a further 2 are newly-detected and spectroscopically 
confirmed in this work. Although three halos have not yet been confirmed, the high success rate within the 
XMM-LSS field (12/15) confirms that weak lensing provides a reliable method for constructing cluster catalogs, 
irrespective of the nature of the constituent galaxies or the intracluster medium. 

Subject headings: cosmology: observations — dark matter — large scale structure of universe 



1. INTRODUCTION 

Clusters of galaxies represent the most massive bound sys- 
tems in the cosmos. Although they result from non-linear struc- 
ture evolution, the departure from linear growth is modest com- 
pared to that for less massive objects. As a result, simple an- 
alytic models can provide an accurate indication of their ex- 
pected number density at various redshifts. This is the primary 
reason why cluster of galaxies are considered to be valuable 
cosmological probes. 

Cosmological attention has focused on the redshift- 
dependent number of clusters, N{z), whose mass exceeds a cer- 
tain threshold. This is one of the most straightforward observ- 
ables, and is a function of the cluster mass function, cfn j dMdz, 
and the evolution of comoving volume dV /dil{z)- The mass 
function is obtained fr om the growth rate of density fluctua- 
tions, S(z), numerically (iJenkins et aLlbOO P) under the assump- 
tion of a particular theory of structure formation, e.g. the 
currently-popular cold dark matter (CDM) model. Since both 
S(z) and dV /dil{z) are dependent on the cosmological model, 
useful constraints could be estimated by comparing N(z) with 
various model predictions. To make progress, e.g. on the dark 
energy equation of state parameter w, data on several thousand 
clusters to z :::^1 i s thought to be required (iLevine et al ] l200l 
IWang et al] |2004|) . and maintaining an accurate and uniform 
mass threshold is critical. 

Most early work focused on selecting clusters optically, with 
detection techniques th at have impro v ed ov er the decades: 
e.g. matched-filter (iPostman et al.l 1 19961). red-sequenc e 
dGladders & Yed 12000 '), cut-and-enhance ( Goto etall l2002h . 
Optically selected samples have traditionally suffered from un- 



certainties in the optical richness - mass relation, although 
there has been recent progress in calibrating the closely related 
richness-velocity dispersion relation using the large maxBCG 
sample of clusters identified in the Sloan Digital Sky Survey 
dBecker et alj|2007l) . 

So far. X-ray sa mple s have been the most popular cosmo - 
logical probes e.g. [g ohringer et all (1200 lb . Ilkebe et all (|2002|) . 
iReiprich & Bohringeri (i2002<) . Luminosity (Lx) or tempera- 
ture (7x)-limited samples offer simpler selection functions be- 
cause the observables, Lx, 7x, are a fair estimate of the cluster 
mass, calibrated through empirical scaling relations. The de- 
rived mass does depend, however, on the assumed dynamical 
state of the system. Unrelaxed clusters, arising for example 
from recent mergers, will introduce scatter in the scaling rela- 
tion. A recent study by Smi th et aTl (l2005l) points out that at 
least half of 10 z 0.2 cluster cores show unrelaxed features 
and a scatter of a ^ 0.4 around the mean scaling relation. 

Weak gravitational lensing, which analyses the coherent 
shear pattern of background galaxies, can potentially provide 
estimates of the cluster mass regardless of its dynamical state 
or the properties of the constituent galaxies. For some years, 
the method has been used to caUbrate m ass obtained from X- 
ray data (Small et al.1ll997l) . lAllen et all (120031) concluded that 
X-ray mass measurements are consistent with those from weak 
lensing, particularly for relaxed systems w hich thus offer a use- 
ful cosmological probe( ^len et alj[2004l) . 

A natural extension of this progress is thus to consider se- 
lecting clusters directly from weak lensing signals. The de- 
velopment of panoram ic imaging surveys ha s now made this a 
practical proposition. IWittman et al.l (1200 lb reported the first 
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discovery of a cluster located from a weak lensing convergence 
map, during the cou rs e of conducting the De ep Lens Survey 
jWittman et alj|2002h . 'Mivaza ki et all (l2002al ) later undertook 
a systematic search of mass concentrations on a 2 deg^ field us- 
ing the Suprime-Cam imager on Subaru. They detected several 
significant (S/N> 4) candidates, one of whic h was later spectro- 
scopically identified as a cluster at z =0.41 . IMetterscheidt et alj 
020051) investigated 50 randomly-selected VLT FORSl fields, 
spanning 0.64 deg^ in total, and found 5 shear-selected candi- 
dates, each associated with an overdensity in luminosity. The 
first results from the Deep Lens S urvey, based on an ar ea of 8.6 
deg^ have also recently emerged dWittman et al.ll2006[) . 

The above pioneering studies have demonstrated that clusters 
can be located directly via weak lensing. However, key issues, 
including the optimum selection threshold, the rate of spurious 
detection and the degree of mass completeness at a given red- 
shift, crucial for any eventual cosmological application, remain 
unresolved. 

At the present time, theoretical studies offer the only insight 
into these issues. Projection is one of the most troublesome as- 
pects of a weak lensing survey, because of the relatively broad 
window function. Unrelated structures contributing to the sig- 
nal would lead to an overestimate of the cluster mass. More- 
over, as the noise in the convergence map arises largely from 
shot noise in the ellipticity distribution of background galax- 
ies, some fraction of genuine clusters might be missed in a 
shear selected cata l og. N-body and ra y -tracing simulations 
jWhite et all 120021: iHamana et all 12004 iHennawi & Spergell 
I2005h have concluded that, for systems whose convergence sig- 
nal lie above a 4 standard deviation (cr) threshold, 60-75% of 
clusters can reliably recovered (completeness). Likewise, for 
peaks detected in the simulated data using typically-used al- 
gorithms, 60-75 % represent genuine clusters (efficiency). The 
difference in these figures between the various studies is largely 
due to differences in the lower mass limit adopted in the studies. 

This series of papers is motivated by the need to address these 
key issues observationally. The survey we desc ribe is a natural 
and ong oing extension of the 2 deg^ survey of iMivazaki et al.l 
(l2002ah : to date a total field of 21.82 deg^ has been imagec . 
This first paper describes the imaging survey and discusses the 
validation of the candidates found, both via spectroscopic ver- 
ification and comparison with clusters located via X-ray tech- 
niques. Later papers in the series will extend the spectroscopic 
follow-up to the full sample and will consider the feasibility of 
deriving cosmological constraints from both this survey and fu- 
ture enhanced versions. W e note that a similar y motivated pro- 
gram h as b een initiated by iMaturi et all (l2006l) . ISchirmer et al.l 
( 120061) and lDietrich eTaiJ ( l2007l) . all of which made used of the 
imaging data taken by 2.5 m VLT survey telescope. We com- 
pare the their conclusions with our own in this paper. 

We n ote that a similarly moti vated program has bee n initi- 
ated by iMaturi et al] (l2006l) and ISchirmer et"aLl (|2006|) whose 
conclusions we compare with our own in this paper 

A plan of the paper follows. We discuss the imaging strategy 
and data analysis in §2, including construction of the cluster 
catalog and its reliability. In §3, we describe our initial spec- 
troscopic follow-up with Subaru and address the completeness 
by comparing X-ray selected clusters on one of our survey area 
where relevant X-ray data is available. We summarize our con- 
clusions in §4. 

2. IMAGING OBSERVATIONS & DATA ANALYSIS 



2. 1 . Survey Fields 

In order to evaluate the efficiency of our weak lensing sur- 
vey for locating cluster halos, we considered that a comparison 
with a samp le of X-ray se lected clusters would be highly ad- 
vantageous (lHenrvll2000h . Therefore, our survey fields were 
primarily selected to contain X-ray data as shown in Table [T] 

We set a minimum ROS AT exposure time, Tr, of 25 ksec en- 
surmg a detection hmit of (0.5-2. 0/teV) - 2 x lO'^^erg/s at 
z ~0.5 (for Hq = 15,^M = 0.3, n\ = 0.7). This corresponds to 
M ~ 10''*Mq, which is well matched to the likely mass detec - 
tion limit of our weak lensing survey (Miv azaki et alj|2b02al) . 
More sensitive X-ray missions, XMM-Newton and Chandra, 
have been surveying ROSAT fields, in part, to deeper lim- 
its. These include the "Lynx", "Lockman Hole" and "UK 13 
hr deep field" in our target list. XMM is also actively in- 
volved in international campaigns of multi-wavelength wide 
field (> Ideg^) observations such as the "COSMOS" , "XMM- 
LSS" and "SXDS" fields which we also inc luded. Among 
these, the XMM-LSS field ( Pierre et al ."2004) is of particular 
interest given its panora mic area and publish ed cluster catalo gs 
(IValtchanov et alJlIool fWillis et al. 2005; P ierre et alJl2006^ ■ 

Finally, we inc luded the DEIMOS DEEP2 survey fields 
dDavis et alj l2003h where spectra of ~ 50,000 faint galaxies 
will eventually become available. This will enable close cor- 
relations of lensing mass and various measures of the luminos- 
ity density as was r ecently pioneered for the COSMOS field 
dMassev et alj|2007l) . Although X-ray data is not available for 
these DEEP fields at the current time, such observations are 
planned and will likely become available soon. 

2.2. Observations 
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Fig. 1. — Seeing (a) and raw stellai' ellpticities (b) measured during the 
observing run of a intensive program (May 1-4 and September 27, 2003). Av- 
erage values derived from 700~1000 stars over the entire field are plotted. 
The abscissa shows the observing time (in days) with the portion of each day 
shifted arbitrarily for clarity. The median seeing is 0.65 it 0.14 arcsec, which 
is typical for the survey. 

The imaging observations were largely carried out on May 
1-4 and September 27, 2003 as part of an "Intensive Program" 
of Subaru Telescope. The Lynx field was observed on January 
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Table 1 

Suprime-Cam Weak Lensing Survey Fields. 



Field 


RA 


DEC 


Area" 


Secure Area*" 


Seeing'^ 


Pgal 














deg2 


deg2 


arcsec 


arcmin^ 


ksec 


ksec 


ksec 


DEEP02 


02:30 


00 


1.39 


0.73 


0.70 ±0.06 


33.5±6.1 








SXDS 


02:18 


-05 


1.12 


0.83 


0.68 ±0.06 


47.7 ±5.7 






100 


XMM-LSS 


02:26 


-04 


2.80 


2.24 


0.55 ±0.07 


46.0 ±6.7 






10 


Lynx 


08:49 


H-45 


1.76 


1.28 


0.80 ±0.08 


30.7 ±7.3 


64 


300 


140 


COSMOS 


10:02 


H-Ol 


1.92 


1 /II 
1.41 


0.54 ±0.03 


in 1 1 '^ 1 

37.1 ±2.1 






30 


Lockman Hole 


10:52 


H-57 


1.85 


1.57 


0.60±0.14 


39.3±7.8 


200 


300 


100 


GD140 


11:36 


H-30 


1.83 


1.50 


0.71 ±0.17 


29.3 ±12.9 


33 






PGl 159-035 


12:04 


-04 


1.43 


1.19 


0.75 ±0.05 


23.4±3.6 


51 






13hrField 


13:34 


H-38 


2.06 


1.72 


0.74±0.17 


29.6±9.6 


110 


120 


130 


GT02deg2 


16:04 


H-43 


2.01 


1.53 


0.67 ±0.04 


38.0±3.6 


26 






CM DRA 


16:34 


H-57 


1.38 


0.99 


0.72±0.12 


28.4±8.4 


47 






DEEP16 


16:52 


H-36 


1.20 


0.93 


0.76 ±0.08 


26.4±4.0 








DEEP23 


23:30 


00 


1.07 


0.80 


0.58 ±0.01 


36.3 ±1.3 








Total 






21.82 


16.72 













"Field area covered by .several Suprime-Cam pointings. 
''Secure area used for halo sample (see text) 
''•''Seeing (FWHM) and galaxy densities pgai refer to average values of the constituent pointings. The scatter is listed in terms of a standard deviation. 
Tr, Tc, and Tn represent X-ray exposure times of ROSAT, Chandra and XMM-Newton, respectively. 



29-30, 2003, SXDS data was obtained from the public archive 
and the COSMOS field observed on February 18 and 21, 2004. 

The Suprime-Cam field size is 0.25 deg^ and all observa- 
tions were conducted with the Rc filter (except the COSMOS 
field which was observed in /'-band to enable direct compari- 
son with ACS /HST's F8 14W images). The total exposure time 
was 30 minutes for each pointing, taken via four 7.5 minute 
exposures in a dither pattern of spacing ~ 1 arcmin. 

Fig. [T] (top) shows the seeing (FWHM) of each 7.5 minute 
exposure over a typical five night observing sequence. The bot- 
tom panel shows ellipticities of moderately bright unsaturated 
stars which provide a measure of the PSF anisotropy; the raw 
ellipticity is mostly 2 - 4 % (and occasionally > 5 %). We dis- 
cuss the question of the PSF anisotropy in Appendix |A] Thus 
far, we have surveyed over 21.82 deg- as summarized in Ta- 
ble [1] 

The useful field area excludes the surroundings of bright stars 
and galaxies and field boundaries. Individual pointings whose 
seeing was worse than 0.9 arcsec were also excluded; this oc- 
curred for only 5 % of the clear time (see section l23T l. 

2.3. Data Reduction 

The data reduction procedures for the presen t survey closely 
followed those described in Section 3.1 of iMivazaki et al.l 
(|2002b), enhanced as discussed below. 

Normally, with Suprime-Cam images, each CCD exposure is 
"mosaic-stacked" to yield a single image of a particular point- 
ing. A simple geometrical model is used for the focal plane 
astrometry whose parameters include the effects of camera dis- 
tortion, the displacement and rotation of each detector from a 
defined fiducial location and the offset and the rotation of the 
dithered exposures. The best fit parameters are obtained by 



minimizing the positional difference of control stars (70^100 
stars per CCD) held common for each exposure. The residual 
alignment error in this procedure is ~ 0.5 pixel rms (0. 1 arcsec 
rms). 

Such a residual is sufficiently small for most imaging appli- 
cations. However, in seeing better than 0.6 arcsec (FWHM), 
a misalignment of 0.1 arcsec between two images introduces a 
~2 % ellipticity on the stacked image which is a serious issue 
for weak lensing studies. 

A further improvement is thus necessary. The residual 
(Ax, Ay) from a reference frame is parameterized as a poly- 
nomial function of field position £= (x,y) as: 

6 / 6 1 

/=0 m=0 /={) m=0 

Pixel values are estimated via linear interpolation of neigh- 
boring four pixels. The coefficients aim,e and bim,e are then ob- 
tained by minimizing the variance of the residuals. A sixth or- 
der polynomial is usually sufficient for this purpose. Each indi- 
vidual CCD image is 'warped' using this polynomial correction 
prior to stacking. This process reduces the alignment error to 
^ 0.07 pixel (0.014 arcsec) and is similar to the "Jelly CCD" 
model described in Kaiser et al. 

Although an external stellar catalog would ideally be used 
for accurate astrometry, no such data is available at the relevant 
faint limits (R > 22). We therefore employ the first exposure, 
coiTected by the simple geometrical model discussed above, as 
the basic reference frame. 

We noticed that this refined procedure still introduces some 
artificial deformation. Fig.|3a) represents the raw image of one 
CCD whereas Fig. |2l;b) is that slighdy rotated by 7 x 10""^ rad, 
a typical value, using the mapping described above. Clearly 
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some degree of artificial deformation is introduced. After some 
experimentation, we found this deformation arises from the un- 
dersampled nature of the Suprime-Cam images. By adopting 
2x2 oversampling prior to rotation, the resuhing ellipticity field 
shows no sign of image deformation. However, this is a com- 
putationally a very time-consuming solution. 
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Finally, we masked all objects close to bright stars (within 18 
arcsec for busNO-A < 15, 90 arcsec for busNO-A < 11.7). Light 
halos around bright stars can introduce spurious galaxies. 
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Fig. 3. — Cumulative galaxy number counts used in the weak lensing anal- 
ysis. Three representative cases in the GD140 field are shown to demonstrate 
how seeing affects the surface density. 
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Fig. 2. — Effect of the under-sampled warp correction in an image taken 
with 0.7 arcsec seeing: (a) no operation, (b) after rotation by 6.8 X 10~* rad. 
A 1.8% residual is reduced to 0.75 %. 

Accordingly, in our final analysis, instead of using oversam- 
pled images, we modified the mapping procedure itself. We 
estimate the pixel values of the target images from 3rd order 
bi-linear polynomial interpolation of 4x4 source pixels|f] rather 
than the linear interpolation four neighboring pixels in the pre- 
vious procedure. This mapping process avoids introducing im- 
age deformation and is significantly quicker computationally. 

2.4. Galaxy Catalogs 

Object finding and shape measurement was executed on the 
mosaic-stacked images using the imcat software suite devel- 
oped by Nick Kaiser A threshold mm=10 was adopted. Photo - 
metric calibration used L andolt standard stars (Landolt II 19921) 
and the faint standards of ( iMaiewski et al.|[T994 ). We adopt the 
Vega magnitude system in the following. 

Galaxies are distinguished from stars via their half light ra- 
dius, r/i, viz: 



where and <Jr* are the half light radius of a stellar image 
and its rms respectively. The galaxy size distribution is shown 
in Fig. H 

Fig. [3] shows the cumulative number density of galaxies as 
a function of Rc-hand magnitude. The surface density exceeds 
50 arcmin"^ when the seeing is superb (0.47 arcsec) and is ~ 15 
arcmin"^ in those poor seeing images (>0.9 arcsec) discarded 
from our analysis (Fig.|3a)). Table [T] lists the seeing and the 
galaxy density for each field. 

^ In the actual implementation, we employ the Numerical Recipes code polin 2 )Pressetal.ll99l . 




Fig. 4. — The size distribution of faint galaxies for the magnitude range 
23 < Rc < 26 observed under under three seeing condition. The size is es- 
timated by circular gaussian FWHM here. Those whose size are larger than 
seeing size are adopted in the galaxy catalogs. 
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Fig. 5. — (a) Seeing dependence of faint galaxy (23 <Rc < 26) surface den- 
sity for all survey pointings except those in the COSMOS field, (b) Pre-seeing 
shear polarisability tensor, Pry , averaged over the galaxies on each pointing ver- 
sus the seeing. The seeing dependence is satisfactorily small for the selection 
threshold adopted. 



2.5. Weak Lensing Analysis 

2.5.1 . Shape Measurements 

Object shapes are represented by the elHpticities, e = 
(ei,e2) = {/ii - /22,2/i2}/(/ii + /22) where Iij are Gaussian- 
weighted quadrupole moments of the surface brightness dis- 
tribution. The point spread function (PSF) of the images is 
usually smeared by various instrumental effects such as opti- 
cal aberrations and the tracking error of the telescope. The PSF 
anisotropy is estimated based on images of stars, and the galaxy 
images are corrected so that images of neighboring stars are re- 
circularized. Galaxy ellipticities are then corrected as: 



(2) 



where the asterisk designates a stellar value, Ps„ is the smear 
polarisability tensor and is mostly diagonal (Kaiser et al. 1995). 
(P~',e)* is evaluated using stars in the field of view and modeled 
as 5th order bi-polynomial function of position. Eqn|2]then ap- 
pUes this for the galaxy images. This correction is carried out 
independently on each pointing. We further justify the correc- 
tion procedure in Appendix lAl 

2.5.2. Shear Estimate 

The shear induce d by gravitational lensing , 7, is diluted by 
atmospheric seeing. iLuppino & Kaisei] (1 19971) developed a pre- 
scription to convert the observed ellipticities to a 'pre-seeing 
shear' as 



1 -c 



(3) 



where is the pre-seeing shear polarisability tensor defined 



Psh is the shear polarisability tensor defined in iKaiser et al.l 
(Il995h . and Pj*, is the stellar shear polarisability tensor. 

Note that (the inverse of) P^ represents the degree of dilu- 
tion. Since the f and P^m are mostly diagonal, we replace the 
tensors in Eqn. |4] with their trace and evaluate P.^ as a scalar. 
The average value, < P^ >, over all galaxies (23 < Rc < 26) of 
each pointing is shown in FigEJb). <Pj> decreases slightly as 
the seeing worsens but the change is not very large (0.4 to 0.3). 
This is because we only select larger galaxies compared with 
the seeing size. Thus, the dilution factor is 30^40 % regardless 
of the seeing. In the mean time, we compare the first compo- 
nent of galaxy ellipticities, ei, of Suprime-Cam and ACS/HST 
images taken in the COSMOS field, and the result is shown in 
Fig|6] The ellipticities observed by Suprime-Cam are in fact di- 
luted by 36 % compared with those of ACS, which is consistent 
with FigHb). 
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Fig. 6. — Comparison of the first component of ellipticity, e[ , of galaxies de- 
tected both on HST/ACS (F814W, 34 minutes) and Suprime-Cam (i'-band, 20 
minutes, 0.54 arcsec seeing) images in part of the COSMOS field (10 arcmin^). 
The best fit slope and the error are also shown. 

To calculate Pj as a function of position we employed a 
"smoothing" scheme ( Van Waerbeke et al] 120001 : lErben et akl 
1200 U lHamanaetani2003h . We took the median P^ for 20 
neighboring galaxies on the r^-magnitude plane (where rg is 
a measure of object size adopted in the imcat suite). In deriving 
the mean, the weight w on an individual measure is taken to be: 



(5) 



+ Q!^ ' 



where is the variance of the raw 7 of those 20 neighbors, 
obtained using the raw P~^. a is the variance of all of the galax- 
ies in the catalog (~ 0.4). In general, the weighted value of a 



quantity < A > is calculated as < A >= I]^[W,A,/I]^ 



^7 ~ ^sh 



(4) 



The method we adopt is based on that adopted by 
iKaiser et al.l d 19951) . More sophisticated methods have 
since been developed and the variants are summarized by 
iHevmans et al.l (l2006h . In their notation, our method is very 
similar to the procedure termed "LV". 
Based on the results of the STEP simulation study, 
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iHevmans et alj (l2006l) concluded t hat b oth the "KSB+" 
method, modified by [Hoekstra et al] d 19981) and im plementa- 
tions of "BJ02" method dBernstein & Jarvis 1 12002|) are able 
to reconstruct input shears to a few percent level. To cali- 
brate our method, we analyzed the simulated data provided by 
iHevmans et alJ (l2006l) . For the model designated "PSF3", we 
underestimate the input shear by 5% for 7,„,(, = 0.1 and 0.05, 
whereas for ^,n,e < 0.01, the difference, jobs-Jme is insignifi- 
cant. An error of 5 % in the recove red shear is competitiv e with 
most of the methods discussed by IHevmans et al.l (l2006h 1 
% is a typical error). A 5% shear error would induce a similar 
uncertainty in the mass estimate of a typical halo. Such an error 
is considered adequate for the applications envisaged. 

2.5.3. Kappa Map 

The dimensionless surface mass density, k{6), is estimated 
from the shear field 7(6') by the lKaiser & Squiresl(ll993h inver- 
sion algorithm as: 



K{e) = - I d^e'Re[D*(e-e')jie')], 

TT 



where D(6) is defined as 



-1 



(01-/^2)2' 



(6) 



(7) 



which is a complex convolution kernel for k to obtain the 
shear 7. In the actual implementation, we smoothed to avoid 
the effect of noise. A smoothing scale of 9c, = I arcmin wa s 
chosen following the discussion given bv iHamana et al.l (l2003l) . 

We adopt 15x15 arcsec^ square grids and then calculate k(6) 
on each grid using Eqn|6]to obtain the k map. In order to esti- 
mate the noise of the k field, we randomized the orientations 
of the galaxies in the catalog and created a Knoise map. We 
repeated this randomization 100 times and computed the rms 
value at each grid point where k is computed. Assuming the k 
error distribution is Gaussian, this rms represents the l-a noise 
level, and thus the measured signal divided by the rms gives the 
signal/noise ratio, 1^(0), of the k map at that point. 

The Gaussian-smoothed signal/noise map is then searched 
for mass concentrations. Hennawi & Spergel (2005) concluded 
that a 'truncated' NEW filter applied to the aperture mass map. 
Map, is the most efficient detection technique. Such optimiza- 
tion may be necessary to improve the efficiency of future very 
wide field surveys where thousands of clusters are sought. Here 
we adopt a simpler technique in order to evaluate its effective- 
ness in a direct comparison of lensing and X-ray techniques. 

2.6. Halo Catalog 

Figure |7] shows the results of a halo search in one of our 
fields: the XMM-Wide field. The red contour shows the k S/N 
map where the threshold and increment are set at 2 and 0.5, re- 
spectively. The blue contour shows the surface number density 
of moderately bright galaxies (21 < Rc < 23). Local peaks are 
searched on the k S/N map and their positions are marked as 
open and filled circles. Figure [TOl to 1211 show the k S/N maps 
for the remainder of our survey fields. 

While visually inspecting galaxy concentrations around the 
detected halos, we noticed that less concentrated halos tend 
to occur preferentially near bright stars and field boundaries. 



Since regions near bright stars are masked (section 12.41 ) con- 
ceivably the discontinuity in the faint background galaxy distri- 
bution could cause spurious peaks. To avoid this, we reject ha- 
los occurring within a 4 arcmin radius of bright stars (busNO-A 
< 11) and within 2.3 arcmin of the field boundary. These re- 
strictions reduce the survey area by 23 % to what we will refer 
to as the secure survey area (16.72 deg^). Halos found within 
the secure area are termed the secure sample. It is certainly pos- 
sible that a significant fraction of halos lying in the non-secure 
area are genuine clusters. We will discuss this further in a later 
paper concerned with their spectroscopic follow-up (Green et 
al, in preparation). 

Detailed inspection of the halo candidates and the spectro- 
scopic follow-up discussed below revealed our completeness is 
high to a limiting signal to noise in the convergence map of 
3.69. Table 1221 IZl fists the secure sample with S/N > 3.69. 
In this table, represents the number of moderately bright 
(Rc < 22) galaxies within 2 arcmin, indicative of the galaxy 
concentration. 

Concerning the optimum threshold for the significance, de- 
creasing it will increase the halo sample but likely introduce 
more false detections. The optimum value should be set based 
on the spectroscopically-observed true/false rate. Investigating 
the rate is a major goal of our study. In this work, the least 
significant spectroscopically-identified halos have S/N = 3.69 
(XMM-Wide n=23). We adopted the threshold of 3.69 for this 
work so that all the spectroscopically followed-up samples are 
included in the table. 

We adopted significance maps for the selection of candidates 
rather than kappa maps. This is because we would like to min- 
imize contamination by the false peaks. However, the effective 
kappa threshold varies over the field, so we may encounter a 
"completeness" problem; i.e. halos that have high kappa value 
are lost from the list. We investigated such omissions in the 
XMM-Wide field (Figure |7]i, and only one halo is found in this 
category, with S/N < 3.69 & k > Kthres where the Kthres is cal- 
culated as Kthres = 3.69 X Noise global- The ''Noise global" is es- 
timated globally over the entire XMM-Wide field kappa map, 
and is 0.018 here. This halo is lost because the local noise is as 
high as 0.022. 

In practice, it will be very hard to generate completely uni- 
form data sets over the entire field of a survey because weather 
and seeing conditions will vary. We will have to optimize the 
Kthres on a field by field basis based on the data quality of each 
field. Therefore, our strategy is the following: at first, we col- 
lect reliable cluster samples based on the significance, and then 
users of the catalog can set their own kappa threshold or mass 
threshold to carry out their studies. This work represents the 
results of the first step above. We list the kappa values in Ta- 
ble |2l6l 12^ for reference. 
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Fig. 7. — The K-S/N map for the XMM-wide field (thick red contour). The lowest contour is set at S/N = 2 with the increments of 0.5. The smoothing scale, 
0(3 = 1 arcmin. The thin blue contours show the surface number density of visible (21 < Rc <23) galaxies. Positions of detected halos (S/N >= 3.693) are designated 
by small circles. Solid circles represent secure halos (see text) whereas unfilled circles represent samples detected outside the secure area. Squares show the location 
of X-ray selected clusters of galaxies published in iValtchanov et al.. 1 ,2004. ). .Willis et al., 1 ,2005) and , Pierre et aL 1 . 2006.) . (VWP data). Large open circles show halos 
followed up by FOCAS and LRIS. 
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Table 2 
List of shear selected halos 



Field 


n 


ID 


RA 


DEC 


kS/N 






tOCAS 


Known'^ 


NEDCj 


Note 


DEEP02 


00 




37.32 


0.63 


4.44 


0.099 


16 






1.35 






01 




37.87 


0.51 


4.39 


0.083 


24 












02 




37.38 


0.41 


4.25 


0.120 


19 






0.73 






04 




37.15 


0.73 


4.11 


0.086 


25 






0.10 






05 




37.31 


0.44 


4.07 


0.131 


20 






1.03 






06 




37.72 


0.69 


3.97 


0.081 


18 






0.86 






07 




37.86 


0.57 


3.95 


0.074 


18 






0.92 






08 


SL J0228.4+0030 


37.12 


0.51 


3.93 


0.085 


49 




0.46(P) 




VOCE 46 




09 


SL J0228.2+0033 


37.07 


0.55 


3.84 


0.102 


33 




0.50(P) 


. 


SDSS CE J037.099808+00.540769 


SXDS 


00 




34.29 


-5.59 


5.33 


0.057 


26 


- 




- 






01 




34.38 


-4.86 


4.14 


0.084 


24 












uz 












1 












03 




34.74 


-4.70 


3.90 


0.044 


19 












04 




34.96 


-5.12 


3.86 


0.069 


21 












05 




34.41 


-4.50 


3.77 


0.056 


26 










XMM-Wide 


00 


SL J0221. 7-0345 


35.44 


-3.77 


8.15 


0.156 


72 




0.43 




XLSSC 006 




01 


CT TAT1 C 1 mil 

JU225. 7-0312 


36.43 


-3.21 


5.72 


0.108 


41 


0.14 






T T) TO „ A 1 /I 

LKIS z = U.14 




02 


CT TAn /I A r\AAC\ 

JU224.4-044y 


36.10 


-4.82 


5.06 


0.074 


40 


0.49 










04 




35.34 


-3.50 


4.91 


0.082 


21 












08 


SL J0222. 3-0446 


35.48 


-3.80 


4.33 


0.081 


29 








LRIS z = 0.41 




10 




36.25 


-4.25 


4.20 


0.062 


23 












12 


SL J0224. 5-0414 


36.13 


-4.24 


4.06 


0.057 


70 


0.26 






T T O A O ^ 

LRIS z = 0.26 




15 


O T T A^ C O l\ A A \ 

SL J0225. 3-0441 


36.34 


-4.70 


3.94 


0.091 


34 


0.26 










16 


SL J0228. 1-0450 


37.03 


-4.84 


3.94 


0.072 


31 


0.29 










17 


O T T A^ C l\ ACW 

SL J022o. 5-0401 


36.63 


-4.02 


3.90 


0.079 


37 




0.34 




"VT O O f A 1 /I 

XLSSC 014 




IQ 


*5T 10997 7 04Sn 


JO. 7^ 


-4 85 


3 81 


064 


43 








J^lcllc cL dl. 




20 




35.98 


-3.77 


3.81 


0.048 


20 












21 


SL J0228.4-0425 


37.12 


-4.43 


3.80 


0.055 


49 


_ 


0.43 


_ 


XLSSC 012 




22 


SLJ0225.4-0414 


36.36 


-4.25 


3.72 


0.073 


43 


0.14 


_ 


_ 






23 


SL J0222.8-0416 


35.71 


-4.27 


3.69 


0.049 


52 


0.43,0.19,0.23 








Lynx 


00 




131.91 


44.80 


5.84 


0.121 


20 












01 




132.59 


44.07 


5.01 


0.083 


43 












03 




131.83 


44.86 


4.57 


0.139 


23 












05 




131.77 


44.85 


4.37 


0.110 


13 












07 




132.69 


44.95 


4.15 


0.105 


31 












08 


SLJ0850.5-I-4512 


132.64 


45.20 


4.02 


0.085 


53 


0.19 


0.24(P) 




NSC J085029+451141,LRIS: = 0.19 




09 




131.47 


44.96 


4.02 


0.076 


26 












10 




133.02 


44.14 


4.00 


0.108 


23 












12 




132.37 


44.38 


3.90 


0.081 


36 












13 




132.41 


44.37 


3.90 


0.072 


31 








Part of n= 12 




14 




132.54 


44.07 


3.86 


0.066 


48 












15 




132.81 


44.35 


3.77 


0.077 


39 












16 




132.31 


44.30 


3.75 


0.072 


44 












17 




131.40 


44.94 


3.74 


0.084 


25 






0.15 




COSMOS 


00 


SL J 1000.7-1-0 137 


150.19 


1.63 


6.11 


0.113 


64 


0.22 


0.20(P) 




NSC J100047-H013912 




01 


SLJlOOl.4-1-0159 


150.35 


1.99 


5.64 


0.098 


32 




0.85(P) 




Einoauenov et al. f20061 




02 


SJ J0959.6-H0231 


149.92 


2.52 


4.74 


0.067 


83 




0.73(P) 




Einoeuenov et al. (2006) 




05 




149.65 


1.55 


3.92 


0.078 


47 










07 




150.19 


2.01 


3.88 


0.070 


36 











"Number of moderately bright {Rq < 22) galaxies around the halo within 2 arcmin 
*redshift obtained by EOCAS MOS 

■^cluster redshift found in literatures (mainly from NED). "P" stands for photometric redshift. 
''redshift estimated from grouping of galaxies whose redshifts are listed on NED. 
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Table 3 

List of shear selected halos (continued from Table [2?61 ) 



Field 


n 


ID 


RA 


DEC 


kS/N 


K 




LUCAS 


Known'^ 


NhDG Note 


Lockman 


00 


SLJ1057. 5+5759 


164.39 


58.00 


6.28 


0.109 


68 


0.60 








03 


SL J1051.5+5646 


162.88 


56.77 


4.97 


0.082 


31 


0.33,0.35 


_ 


_ 




05 


SL J1047.3+5700 


161.84 


57.01 


4.56 


0.103 


56 


0.30, 0.24 


- 


- 




06 


SL J1049.4+5655 


162.35 


56.93 


4.51 


0.095 


47 


0.42 


- 


LRISz = 0.31 




09 


SL J1055.4+5723 


163.86 


57.38 


4.22 


0.086 


20 


- 


- 


LRIS z = 0.38 




10 


SL J1051. 6+5647 


162.92 


56.78 


4.20 


0.068 


49 


0.33, 0.25 




0.05 part of SL J1051.5+5646 




^ J 


jIj J i LfJ J .^+ J / ZU 




J / . 


4 08 


064 


50 




34 


T1 nS'^ ^-i-S71 Q 




12 






J 1 .J J 


4 07 


Of^S 


26 










1 




162 91 


Jo.uz 


4 04 


OSQ 


21 






U.Uo 




14 




162 54 


57 28 




070 


JO 










15 


SL J1048. 1+5730 


162.04 


57.51 


3.89 


0.071 


35 


0.32 








16 




163.85 


57.95 


3.83 


0.072 


24 




- 


0.02 




18 


- 


163.16 


57.88 


3.77 


0.069 


26 


- 


- 


- 




19 


- 


164.21 


57.70 


3.77 


0.070 


22 


- 


- 


- 




20 


- 


163.23 


57.84 


3.73 


0.102 


29 


- 


- 


- 




21 


- 


163.14 


57.82 


3.72 


0.111 


19 


- 


- 


- 


GDI 40 


00 


SLJl 135.6+3009 


173.91 


30.16 


4.98 


0.126 


35 


0.21 








01 




173.89 


30.21 


4.19 


0.086 


23 










02 




173.96 


29.81 


4.09 


0.069 


17 










03 


SL Jl 136.3+2915 


174.09 


29.26 


4.03 


0.100 


24 






LRIS z = 0.20 




05 




174.77 


29.89 


3.86 


0.111 


26 










06 


- 


174.86 


30.33 


3.83 


0.080 


21 


- 


- 


- 


PGl 159-035 


05 


SLJ1201. 7-0331 


180.44 


-3.53 


4.42 


0.077 


49 


0.52 








06 




180.99 


-3.09 


3.90 


0.119 


21 




- 


0.09 




08 


- 


181.76 


-3.27 


3.71 


0.102 


37 


- 


- 


- 


13 hr Field 


00 


SL J1334.3+3728 


203.60 


37.47 


4.33 


0.128 


74 


0.30 


0.48(P) 


NSCS 11 33424+372S22 




01 


SL J1335.7+3731 


203.94 


37.53 


4.10 


0.091 


65 


0.41 








04 


SLJ1337.7+3800 


204.43 


38.01 


3.85 


0.080 


34 


0.18 


- 


- 




06 


- 


203.85 


37.90 


3.78 


0.068 


27 


- 


- 


- 




07 


- 


204.22 


37.54 


3.77 


0.078 


30 


- 


- 


- 


V. 1 1 V / zLie^ 


00 


T1 609 Sj-d'^'^S 

Jl^ J iOUZ.O+^JJj 


240 72 


43 59 


6 65 


110 


56 


42 








n 1 




9 AO 7S 


49 1^ 
tz. / 


S Al 


n 1 06 


S7 






T T?T^ V — 1 H 




uz 




9zL1 S9 
Z+i .OZ 


A^ 1 Q 


S 1 7 
J. 1 / 


1 


9 1 
Zl 










04 




241.95 


43.60 


4.49 


0.093 


17 










07 


SLJl 604. 1+4239 


241.04 


42.65 


4.19 


0.083 


39 






LRIS z = 0.30 




08 




241.63 


43.61 


4.16 


0.093 


25 










09 


SLJl 605. 4+4244 


241.36 


42.74 


4.09 


0.064 


36 


0.22 








10 




241.18 


43.46 


3.85 


0.092 


29 










11 




241.70 


43.64 


3.84 


0.077 


33 










12 


SLJl 603. 1+4243 


240.78 


42.72 


3.82 


0.091 


37 






Partof SLJl 603. 1+4245 


CM DRA 


04 


SLJ1639.9+5708 


249.98 


57.15 


4.15 


0.073 


32 


0.20 








06 


SLJl 634. 1+5639 


248.55 


56.66 


3.97 


0.104 


31 


0.24 






DEEP 16 


00 


SLJ1647.7+3455 


251.94 


34.93 


4.30 


0.084 


42 


0.26 








01 




253.54 


34.98 


4.24 


0.078 


24 










02 




252.36 


35.02 


3.75 


0.123 


25 










03 




251.79 


35.04 


3.72 


0.087 


21 








DEEP23 


00 


SL J2326.4+0012 


351.61 


0.20 


4.41 


0.076 


19 


0.28 








01 




351.75 


0.00 


4.37 


0.095 


20 










02 




352.33 


0.15 


4.27 


0.107 


33 






1.38 




03 




352.47 


-0.06 


4.11 


0.079 


34 






0.07 




04 




352.22 


0.09 


3.87 


0.084 


21 






1.37 




05 




353.09 


0.04 


3.70 


0.066 


22 
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V 



-a 




Galaxy density [arcmin ^] 

Fig. 8. — Faint galaxy density {p^al) dependence of number density 
of secured halos (S/N > 3.69) in each field. The error bai's are based on 
error estimates. Solid lines shows the predicted halo number density. 

The galaxy distribution is assumed ^ 

: 2.0 and (3 -- 



exp 



.r[(l+a)//3)J 

where a = 2.0 and /3 = 1.5 is adopted. The mean redshift, (z) is related to z, 
by (z) =z*r[(2 + a)//3]/r[(l+a)//3] Three different (z> are assumed; 0.8, 
1.0 and 1.2 (from bottom to the top solid line). 

Here we compare the detected halo number density 
[deg'^l with the pre diction of numerical simulations done by 
iHamana et all (12004'). who have attempted to reproduce a sur- 
vey such as ours. Figure[8]shows the halo density of our thirteen 
survey fields. The horizontal axis shows the density of faint 
galaxies {pgai) used for weak lensing analysis. The error bars 
are based only on y/N estimates and do not include the effects 
of cosmic variance. Solid lines in Figure [8] show the prediction 
of the simulation where three different mean redshifts for the 
background galaxies are assumed (from bottom to top: 0.8, 1.0 
and 1.2). 

Although the scatter is large, there is reasonable agreement 
between the observations and the prediction. We also see a 
gradual increase of the observed number density over the pgai 
range sampled as expected. These comparisons validate our 
observational procedures. We have, however, two outliers in 
Figure [8] These could be caused by cosmic variance or may 
arise from some other reasons. Further follow-up studies will 
be important to clarify the issue. 

3. CLUSTER IDENTIFICATION AND VERIFICATION 

Armed with the shear-selected halo catalog, we now discuss 
the tests we have made to verify its reliability, using both spec- 
troscopic observations and comparisons with X-ray data in the 
fields where the overlap of targets can be studied. 

3.1. Spectroscopic Follow-up 

Because the lensing kernel (or window function) has a rel- 
atively broad redshift range, the superposition of multiple low 
mass halos with different redshifts could yield highly signifi- 
cant weak lensing signals. Such 'superposition halos' are un- 
welcome in a mass-selected cluster catalog. In order to identify 
such superposition halos we undertook a more comprehensive 



spectroscopic survey using a multi-object spectrograph for se- 
lected candidates. This also gives us the velocity dispersion 
of member galaxies, which is an estimate of dynamical mass 
of clusters. By comparing dynamical mass with weak lensing 
mass we will be able to discuss the dynamical state of the clus- 
ters (Hamana et al. in preparation). 

We used the FOCAS spectrograph on Subaru whose multi- 
object mode permits the simultaneous observation of 25 - 30 
galaxies in the field of a particular halo. We used the 150/mm 
grating and a S Y47 order sorting filter. This configuration spans 
the wavelength range 4700 -9400 A (' Kashikawa et al.ll2002h . 
Target selection was based primarily on apparent magnitude 
and the color information was used when available. The ex- 
posure time was 45 ~ 70 minutes depending on the magnitude 
of the selected galaxies and the observing conditions. The spec- 
troscopic data was reduced using standard IRAF procedures 
(Hamana et al. in preparation). Figure |9] shows a typical FO- 
CAS observation where the target is identified as a z = 0.6 clus- 
ter. 

Since May 2004 we have observed 26 halos from the secure 
sample with FOCAS. Higher priority was given to more signifi- 
cant halos, except in the early stages of program (for example in 
the Lynx and PG 1159-035 fields) when the follow-up strategy 
was still being evaluated. Each of the 26 halos has been reliably 
identified with a cluster of galaxies. The redshift so determined 
is shown in the column labeled "FOCAS" of Table l2.6l2.6l 

In parallel with the FOCAS follow-up, long slit observations 
with Keck/LRIS have been carried out to enlarge the sample 
of redshifts as quickly as possible. We have to keep in mind 
that this method cannot discriminate the projected samples, and 
some statistical consideration is necessary in dealing with the 
data for further studies. A complete discussion of that aspect 
of our program is discussed in a separate paper (Green et al, in 
prep.). In the "Note" column of Table IZ6I2. 61 we list the pre- 
liminary Keck/LRIS results as "LRIS z = zvalue" whose iden- 
tification is already robust (e.g. at least two galaxy redshifts 
agree) including one of the halos in the XMM-Wide field (SL 
J0222. 3-0446, z=0.41) where we discuss the reliability of our 
catalog (see section l33T l. 

3.2. Correlation with Published Data 

In addition to undertaking our own spectroscopic obser- 
vations, we also searched the NASA Extragalactic Database 
(NED) for clusters and groups with published redshifts close 
to the position of our detected halos. We included recently- 
publi shed cluste r cata lo gs in the XMM- wide survey by 



Valtc hanov et all (|2004 . Iwilhs et al.1 (l2005h and iPierre et all 
(12006) (VWP data hereafter). In making assignments, we con- 
sidered a cluster to be associated with a halo if the angular sep- 
aration was less than 2 arcmin and the 3 dimensional distance 
less than 1 /i"' Mpc. 13 out of 100 halos can be identified in this 
way with rich clusters in the literature and the published red- 
shift is shown in the column labeled "Known" of Table |2.6|2.6l 
Three of these thirteen clusters were also observed by FOCAS 
in multi-object mode and we note that the mean FOCAS red- 
shift is different from the published value for these systems 
(& = 0.02 ~ 0.18). These discrepancies can be explained by 
the fact that NED redshift of these clusters are all photomet- 
ric, which have larger uncertainties. Therefore, we adopted the 
FOCAS redshifts for these three cases. We searched the NED 
galaxy group catalog in the same manner, but no group matched 
any of our halos. 
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Fig. 9. — (a) K S/N contour map of the most significant lialo in tlie Loclcman field (SL J1057. 5+5759) superimposed on the Rr-band image taken by Suprime-Cam. 
Contours start at a S/N = 22 with an interval of 1. Small circles show the positions of observed galaxies by FOCAS with the redshifts obtained, (b) Redshift 'cone 
diagram' of the observed galaxies. 



Next, we searched for individual galaxies with published red- 
shifts within 3 arcmin of our halos. If suitable data was found, 
we grouped them together using a "friend of friend" algorithm 
with a linking length of 1 Mpc. If we found a clustering of 
more than two members, we calculated the average redshift and 
the mean astrometric position. Using the same proximity crite- 
rion above, we assigned a redshift to 14 further halos. These 
are shown in the column labeled "NEDG". 

In order to estimate the probability of a chance coincidence 
of the matching procedures, we first randomized the position of 
the detected halos in the field, and applied the same matching 
procedure described above to the randomized halo catalog. This 
process suggests that the probability of a chance coincidence is 
roughly 10 % and 25 % for the cluster and galaxy clustering 
searches, respectively. Clearly care should be taken in analyz- 
ing such data. The total number of reliably identified clusters 
are 41 (out of 100 in Table [27612. 6t where we do not include 
those identified by "NEDG" galaxies because the chance coin- 
cidence is not negligible. We assign ID labels only for these 
reliable halos in the tables. 

3.3. Reliability of the Shear-Selected Halo Catalog 

We now turn to the important question of the evaluating the 
reliability of our shear-selected catalog. We will do this by 
examining both the success rate of our identifications and by 
comparing with X-ray samples obtained via the XMM-Wide 
field survey. VWP have confirmed spectroscopic redshifts for 
28 X-ray selected clusters in the 3.5 deg^ survey field. Their 
locations are plotted as squares on Figure [T] Our secure sam- 
ple of 15 shear-selected halos with S/N > 3.69 is plotted as 
solid circles. Those we have spectroscopically confirmed using 
Subaru's FOCAS and/or Keck LRIS are marked by large open 
circles. We find that three halos (1, 16 and 23) are confirmed as 
clusters by FOCAS and are not reported in VWP. 



Among our 15 shear-selected halo samples, only three (4, 10 
and 20) have yet to be identified. This identification success 
rate (80 %) can be regarded as a lower limit given future obser- 
vations may yet locate an associated cluster Although affected 
by small number statistics, this minimum e fficiency is already 
higher than expected by simulation studies( yhite et all 120021: 
iHamana et ani2004t iHennawi & SpergeI.2005lK 

It is interesting to compare our identi fication success rate 
with those found in the GaBoDS survey. iMaturi et al.l (l2006h 
found 14 significant halos in their 18 deg^ area. Among them, 
5 halos turn out to be known clusters of galaxies, 2 seem to 
have associated light concentrations but no spectroscopic con- 
firmation, and the remaining uncertain 7 (50 % ) had no appar- 
ent counterpart in either optical nor X-ray data. Schirmer et al.l 
(2006) undertook a search on the same data adopting a different 
peak selection algorithm and a lower detection threshold. They 
found 158 "possible mass concentration" on the 18 deg^ field. 
If those halos are divided into the same classes as a bove (known, 
light, uncertain), the ratio is almost the same as 'Maturi et al.l 
(i2006 ). Regardless, almost half of the GaBoDS halos could be 
considered uncertain at this point. 

Meanwhile, we find 15 halos in a 2.24 deg^ field (XMM- 
wide) and show that 80 % of have been already identified as 
clusters. It is too early to make any definite conclusion be- 
cause their spectroscopic is underway. However, our (tenta- 
tively) higher success rate could be explained at least in part 
by the larger number density of faint galaxies, pg^i, usable for 
the weak lensing analysis owing to larger aperture and better 
average seeing. In the case of GaBoDS, pgai spans from 6 to 
28 arcmin"^ depending on the field; the average value is 11 
arcmin"^. This is generally smaller than our pgai shown in Ta- 
ble [T] As a result one can expect a lower angular resolution of 
the K map, reducing the S/N ratio for a fixed smoothing scale 
and possibly an increased degree of contamination in the result- 
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ing halo catalog. Another possibility to explain the increased 
contamination is that their sample consists of a combination of 
different sets of catalogs, each of which is selected by different 
methods. This could decrease the significance threshold effec- 
tively, and could introduce more false peaks. 

3.4. Superposition of Multiple Clusters 

iHamana et a D ('2004) estimated, based on their simulation, 
that the halo superposition rate in survey such as ours should 
be roughly 3 % - a small but not negligible effect. A longslit 
spectroscopic survey, such as that undertaken with LRIS (Green 
et al, in prep) might be poorly-suited for locating such cases. 
However, the FOCAS multi-object survey reported here should 
reliably find them. In fact, we have found three apparent su- 
perposition halos (SL J0222.8-0416, SL J1051.5H-5456 and SL 
J1047.3H-5700) out of 26 examined. The overlap rate is broadly 
consistent with expectation considering the small number so far 
sampled. 

4. CONCLUSIONS AND FUTURE PROSPECTS 

We have introduced a new Subaru imaging survey and de- 
scribed techniques for locating and verifying shear-selected ha- 
los. Across a search area of 16.72 deg^ we have found 100 halos 
whose K S/N exceeds 3.69. We have described the first phase 
of a detailed follow-up campaign based on multi-object spec- 
troscopy of 26 halos using FOCAS on the Subaru telescope. 
A later paper in this series (Green et al, in prep) will extend 
the spectroscopic survey to the full sample using a longslit ap- 
proach. 

Detailed studies on one of our fields, the XMM-wide field, 
show that 80% of the shear selected 15 halos in the 2.2 deg^ 
area can be confirmed as genuine clusters of galaxies. 10 over- 
lap with X-ray detections and two are new systems confirmed 
spectroscopically. The overall success rate and reliability of our 
sample provides convincing proof that, with care, a weak lens- 
ing survey can provide a large sample of mass-selected halos. 

We compare our success rate and the reliability of our catalog 
with that of the GaBoDS survey and conclude a major advan- 



tage of our approach is the superior imaging depth which leads 
to a high surface density of usable galaxies. This suggests fu- 
ture, more ambitious, surveys for shear-selected halos will be 
more effective if undertaken with large aperture telescopes. 

It is interesting to use our results to estimate the requirements 
for a future survey motiv ated by the need to constrain dark en- 
ergy. iKolb et aTl (l2006l) discuss hypothetical missions which 
would aim to analyze the redshift distribution N{z) of 10,000 
clusters. Figure |8] shows that our survey technique typically 
finds 5 halos deg"^. Because the field size of Suprime-Cam is 
^ 0.25 deg^, it takes two hours to cover 1 deg^ assuming the 
exposure time of 30 minutes of each field as adopted in this 
study. A survey of 10,000 clusters would be prohibitive even in 
terms of imaging alone (~ 500 clear nights) even before con- 
templating the follo w-up spectroscopy. T he Hyper Suprime- 
Cam (HSC) project (Mivaza ki et al.ll200 6') has been proposed 
to remedy this shortcoming. HSC is expected to have a field 
of view ten times larger than Suprime-Cam while maintaining 
the same image quality. This new facility will makes the 2000 
deg^ scale imaging survey within a reasonable number of clear 
nights. 

Spectroscopic follow-up might be enabled by the proposed 
prime focus multi object optical spectrograph WFMOS (Wide 
Field Fiber Multi Object Spectrograph) whose field is 1.5 deg^ 
in diameter Typically we can expect ~10 shear-selected clus- 
ters in each spectroscopic field. With only 20 targets per halo, 
only a small fraction of the several thousand fibers envisaged 
for WFMOS need be allocated to the halo verification program. 

We are very grateful to Subaru astronomers: Y. Oyama, K. 
Aoki and T. Hattori for their dedicated supports of the FOCAS 
observing. Numerical computations presented in this paper 
were carried out at the Astronomical Data Center (ADC) and 
at the Center for Computational Astrophysics (CfCA) of the 
National Astronomical Observatory of Japan. This work is sup- 
ported in part by Grant-in- Aid for Scientific Research (Kaken- 
hi) of Japan Society for the Promotion of Science (JSPS): 
Project number 15340065 (SM&TH) and 17740116 (TH) . 
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Fig. 10. — DEEP02 field k-S/N map (red contour) and surface number density of moderately bright (21< i?c <23) galaxies (blue contour). Legends are the same 
as Fig|2] 
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Fig. U.— 



SXDS field ft-S/N map and surface density of galaxies, (blue contour). 
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Fig. 12. — Lynx field k-S/N map and surface density of galaxies. 
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Fig. 13. — COSMOS field k-S/N map and surface density of galaxies. 
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Fig. 14. — Lockman field k-S/N map and surface density of galaxies. 




Fig. 15. — GD140 field K-S/N map and surface density of galaxies. 
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Fig. 16. — PG1159-035 field k-S/N map and surface density of galaxies. 
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Fig. 17. — 13 hr field k-S/N map and surface density of galaxies. 
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Fig. 18. — GTO 2 squai'e degree field k-S/N map and surface density of galaxies. 




Fig. 19.— 



CM DRA field k-S/N map and surface density of galaxies. 



18 



Miyazaki et al. 



35.2 

O 

O 35 

u 

Q 34.8 



254 




® a 



353.5 



253 252.5 

RA2000 



Fig. 20. — DEEP16 field k-S/N map and surface density of galaxies. 
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Fig. 21. — DEEP23 field k-S/N map and surface density of galaxies. 
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PROBABILITY TO OBTAIN SPURIOUS PEAKS DUE TO INSUFFICIENT CORRECTION OF THE PSF ANISOTROPY 
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The shear induced by massive cluster of galaxies is expected to be 7~ 10 %. This is actually larger than raw ellipticities due to 
PSF anisotropy (2 ~ 4 % refer to Fig. [TJ, but the difference is not so significant. Therefore, the correction of the anisotropy is very 
important to make a precise kappa map. Here, we examine the PSF anisotropy Suprime-Cam, and evaluate the effect of imperfect 
correction in this work. 




Fig. A22. — Two examples of ellipticities field. Orientation of each bar shows the direction of the major axis, and the length scales as its ellipticity. No subtraction 
of the offset is made. Arrow on left top of the panel shows the size of 10 % ellipticity. The seeing was 0.65 arcsec (FWHM). 

In order to characterize the anisotropy of the Suprime-Cam images, we obtained sequences of i'-band short exposures (1 min) of 
dense stellar fields with various telescope pointings over one night long. The seeing was mostly ^ 0.65 arcsec (FWHM). The PSF 
anisotropy is approximated as an elliptical, and the field position dependence of the ellpticities are investigated. Two examples of 
such ellipticity fields are shown in Fig. lA22l General tendency that we note is that the fields can be represented as a super-position 
of (1) the radial field at four corners (almost invariant) and (2) almost unidirectional field (variable). The variable components is 
most likely due to the telescope shaking whereas the invariant component can be explained by the optical aberration of the corrector 
Slight asymmetry of the radial component is visible (i.e. ellipticities near the lower left corner is larger than other corners), and this 
would be a sign of imperfect optical alignment. 
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Fig. A23. — (a) Ellipticities of stai's azimuthally averaged over the field of view (0.65 arcsec seeing i'-band 60 sec exposure). Open circle and squai'e show the 
results of "before" and "after" the anisotropy correction, respectively (see text). Sold squares show expected ellipticities from perfectly aligned optics under the 
seeing of 0.65 arcsec. (b)Peak distribution function of kappa S/N map created by residual of anisotropy corrections. 

We also note that the discontinuity of elUpticities is not visible across the boundary of the CCDs, and the ellipticities can be 
represented as a single continuous function of field position. We adopt 5-th order polynomial function here. In order to simulate 
the actual science field analysis, we randomly select controls stars and "galaxy role" stars from the star catalogs with an appropriate 
density; 2 arcmin"' and 40 arcmin"\ respectively. Using the control stars, we obtain the best fit coefficients of the polynomial, and 
the anisotropy of "galaxy role" stars is corrected using Eqn. |2] Because the intrinsic ellipticities of "galaxy role" stars is all zero, the 
residual ellipticities after the correction is a estimate of incompleteness of the correction. Fig |A23r a) shows the azimuthally averaged 
ellipticities of stars; before the correction, the typical 3 % raw ellipticity is seen. It might be interesting to note that the ellipticities of 
3 %. can be induced by rms pointing error of merely 0. 1 arcsec under the seeing of 0.65 arcsec. The ellipticities is reduced down to 
about 0.75 % after the correction. Beyond the field angle of 18 arcmin, the correction does not work fine and we decided to eliminate 
the field r > 18 from lensing analysis; which, however, results in only a few percent loss of FOV. 

We estimated a PSF anisotropy caused by the optical aberration using a ray-tracing code, zemax. The calculated PSF is convolved 
with 0.65 arcsec FWHM gaussian, and the shape is evaluated by elliptical. The result is shown in filled square in Fig IA23f a). 
Compared with this ideal case, the observed ellipticities is large even after the correction. This shows a limit of the adoption of 
single polynomial function as a representatives of the PSF anisotropy, where the residual is still locally correlated correlated and not 
completely random. We now want to evaluate the impact of the incompleteness onto the k S/N map. The map is create d based by 
the re sidual ellipticities after the correction. We simulate the galaxy ellipticities using the following "conversion" formula dHoekstra I 
l200l . 

pSal 

e«"' = -^e*, (Al) 

which is essentially a sensitivity correction against the PSF anisotropy. This is necessary because galaxies are larger and their shape 

is more insensitive to the deformation compared with stars. We evaluate -ft with )J|.^p'" where we adopt {TriPlm )) = 0.1 which is 

a typical value under the typical 0.7 arcsec seeing. Since the {Tr(P*^)) ^ 0.2 here, the conversion factor is roughly i. We omit 
correction because it cancels out in this case when we calculate the S/N. We created twelve such maps from independent exposures, 
and co-added the peak distribution functions to obtain Fig |A23r b). It is obvious that the incompleteness of the anisotropy correction 
is quite unlikely to create any significant (say S/N > 3) fake peaks. 



